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(Dated: April 3, 2005) 

We study the driven tunneling of a one-dimensional charged particle confined to a rectangular 
double-well. The numerical simulation of the Schrodinger equation based on the Cranck-Nicholson 
finite-difference scheme, shows that the modulation of the amplitude of the external field may result 
in the parametric resonance. The latter is accompanied by the breakdown of the quasi-periodic 
motion characteristic of the usual driven tunneling, and the emergence of an irregular dynamics. 
We describe the above breakdown with the occupation probability for the ground state of the 
unperturbed system, and make the visualization of the irregular dynamics with the help of Shaw- 
Takens' reconstruction of the state-space. Both approaches agree as to the values of the resonant 
frequency for the parametric excitation. Our results indicate that the shape of the laser pulse could 
be essential for generating chaotic tunneling. 
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1. Driven transitions in a double well are instru- 
mental for studying the tunneling in various fields of 
physics and chemistry, 0- Considerable attention 
has been drawn to the tunneling dynamics in the pres- 
ence of a driving force with a time dependent amplitude. 
In his seminal paper, [3], M. Holthaus showed that shap- 
ing the driving force may be instrumental in controlling 
the tunneling in a bistable potential. In particular, it 
was shown, Q, that choosing an appropriate envelope 
for laser pulse, one may perform the population trans- 
fer on time scales much shorter than the base tunneling 
time. This situation is intimately related to the prob- 
lem of quantum chaos, which is generally approached 
within the framework of the quasi-classical approxima- 
tion and Gutzwiller's theory. In fact, classical chaotic 
systems are often used as a clue to the quantum ones. 
In contrast, it would be very interesting to look at the 
quantum chaos the other way round and consider sys- 
tems which need studying without approximations that 
could have bearing upon classical mechanics, for exam- 
ple, particles confined to potentials of a size comparable 
with the de Broglie wave length. This has also an addi- 
tional interest owing to the fact that calculations within 
the framework of semiclassical theory should depart from 
the quantum ones on the time scale of h/AE, AE being 
the typical spacing between energy levels. Therefore, the 
Schrodinger equation describing the problem needs nu- 



merical studying. 

This letter aims at the theoretical description of the 
tunneling transitions for a charged particle in a one- 
dimensional rectangular double- well potential U (x) and 
an external periodic field V(t) that mimics electro- 
magnetic irradiation. The Hamiltonian of the system in 
dimcnsionlcss units is given by the equation 

H = \tf + U(x) + V(t) 
in which the potential U reads 



—Ui, a < x <b 
U = I -U 2 , c<x <d 
0, otherwise 



(1) 
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i.e. the unsymmetrical double well comprising the two 
shafts of different depths Ui,U%, see FIG^ The shape 
of U is chosen in such a way that the values of the wave 
functions, ipo(x,t) =< x\0 >, ip\{x,t) —< x\l >, of the 
ground and the first excited state, are very small outside 
the right and the left shafts of the well, respectfully, see 
FIGH Consequently, the particle's transitions |0 >^ 
1 > can be visualized as the tunneling ones. The field 
frequency f2 verifies the constraint of resonance, which 
in dimensionless units reads fi = Ei — Eq , where Eq , Ei 
are the energies of the ground and the first excited states. 
The external field V(t) has the form 

V(t) = - A sin(fit)j3, p = -id x 
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FIG. 1: Rectangular one-dimensional double-well potential. 
Lines "0" and "1" indicate the moduli of the wave functions 
|i/>o(£,t)| and |^i(a;,t)| of the ground and the first excited 
states. In the chosen configuration D — 0.876 is the distance 
between the wells; Lq — 2.337, L 1 = 2.045 are the widths, 
and Uo = 13.82, Ui = 11.91 are the depths of the shafts in 
dimensionless units. 



where p is the operator of momentum, corresponding to 
an electro-magnetic wave described by the vector poten- 
tial A x — A z — 0,A y = A sin £l(t — x/c), the scalar po- 
tential cf> being equal to zero, [5j- In what follows we 
neglect the dependence of A on x since the wave length 
is assumed to be much larger than the well's size. 

We are looking for regimes in which the regular tun- 
neling breaks down and becomes chaotic. The main in- 
strument employed for this end, is the Shaw-Takens (ST) 
method, 3, for reconstructing the state-space of a sys- 
tem from a stream of data, obtained with the numerical 
integration of the Schrodinger equation, by embedding 
it in a phase space of enough dimension T>, in our case 
T> = 3. The important thing is that one may monitor the 
dynamics of driven tunneling, by using the occupation 
probability iVo of the ground state of the unperturbed 
system 

TV = | <0|^> | 2 , 

where |0 > and \ip > are the ground vector, and the 
state vector of the system at the moment of time t, re- 
spectively. To obtain the visualization we shall form a 
time series that comprises the values of No at moments of 
time nr, n = 0, 1, 2, . . . ,; r being the lag-time, so that 
Nok = No(kr), k = 1, 2, 3, . . . . Then the series of vec- 
tors Yfc = {N 0k ,N 0k+1 ,N Qk+2 ),k = 0, 1, 2, 3, . . ., serves a 
3d- visualization of the state-space for the given problem. 
In particular, we can observe the transition from the reg- 
ular motion characteristic of the unperturbed resonant 
motion to the irregular one generated by a parametric 
excitation given by Eq.@ below, see FIG[21 

2. The parametric excitation of the system is gener- 
ated by the modulation of the amplitude A of the vector 



A = A [l- esin(wt)] 



(2) 



The approximate expression for the frequency of the 
parametric resonance caused by can be obtained by 
employing the rotating wave approximation, |fj after 
some analytical calculations. It reads 



u prm 



= Ak, 



< 0\p\l > 



(3) 



In fact, the actual resonant frequency uj prm differs from 
that given above and its calculation requires some nu- 
merical work. 

3. The visualization with the help of the ST-method 
gives a torus in 3d-space, fairly well drawn, changes in 
the lag-time resulting only in continuous deformations of 
the picture, and no sudden modifications observed, see 
FIG|21 In contrast, in case the parametric excitation at 
a frequency close to the value given by being present, 
the choice of the lag-time is very important for obtaining 
a meaningful visualization. The break down of the quasi- 
periodic motion by the resonant parametric excitation is 
shown in FIGJ2J It is to be noted that the visualization 
depends on a number of premises, and first of all the value 
of the lag-time r and the dimension d of the visualization 
window. The wise choice of r is dictated by characteristic 
times of the motion under the investigation, so that for 
certain values of r the visualization picture is coherent 
enough whereas for others it is not meaningful. There- 
fore, using the ST-method one has to compare pictures 
obtained for different values of the lag-time. 

4. We may look at our problem the other way 
round. Consider sets of the system's states Fkk+ii & 
1.2 TV, where TV is large enough, given by constraints 



A k < N < A k+1 , 



(4) 



where are intervals dividing segment < A^o < 1 in 
Af equal parts. Take a period of time T, sufficiently large, 
and consider times r k spent by the system in the sets of 
states .Ffcfe+i, that is where (0J is verified. We shall define 
the visiting frequencies as 



T 



One may cast the v kk +\ in the form of the probability 
density, £(iV ), for the distribution of iVo, by assuming 
that all the intervals A^ +1 — A k be of equal size A, and 
defining 



A 



On considering the limit of £ as TV — > oo, we shall get the 
probability density f (JV ), see FIGEJ 

To put it in a quantitative analytical form, we may 
consider the characteristic function Xab(x) = 1, if a < 
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x < b and otherwise, and introduce the quantity 

Ma, b) = Yrn^ ji J Xab(No) dtj , (5) 

so that Vab can be considered as the frequency of visiting 
a region of states determined by constraint a < < b. 

The important point is that the limit indicated in the 
above equation does exist in the context of the prob- 
lem under investigation. To see the fact we shall employ 
the normalization of wave functions in a finite box, that 
is, our system being confined to a rectangular potential 
of infinite depth of a size much larger than the double 
well's size. Then the system has the discrete spectrum 
of eigenvalues, and the characteristic function Xab{No) 
can be expanded in Fourier series, the time averaging in 
the above equation results in cancelling out terms oscil- 
lating in time, so that the limit exists. The use of the 
finite box normalization is also important for managing 
the numerical simulation. It is alleged to be known, the 
finite integration mesh, or basis, results in reflection of 
the wave packet that distorts the time evolution of the 
packet. By imposing the finite box normalization we use 
the physical framework that appears to be compatible 
with the tunneling dynamics in the double well of a size 
much smaller than the box and enables us to overcome 
the numerical artifacts. 

It is worth noting that the construction of visiting fre- 
quencies is analogous to that of the dwell time which is 
aimed at studying the localization of wave packets and 
defined by the equation, 0|, 0, 

/+oo r-b 
dt I \i/j(x,t)\ 2 dx 
-oo J a 

The probability density £ indicates that at the reso- 
nant frequency u prm the character of the system's motion 
undergoes a drastic change, see FIG|3 At this point it 
should be noted that the frequency of the parametric res- 
onance Lo prm which can be found by comparing the sizes 
of the peaks in FIG|21for different values of u> in equation 
, differs by approximately 5% from the value given by 
the rotation wave approximation @ . The resonant value 
corresponds to the most pronounced breakdown of the 
twin peaks, which correspond to the ground and the first 
excited state of the free system. 

5. Concluding we should like to note that, as follows 
from Eq.JSJ) for the amplitude of the driving force, the 
pulse V has a triplet structure determined by the main 
contribution at frequency Q = E\ — Eq and two satellites 
at fi ± Up rm , of less amplitude, owing to the equation for 
A(t) given by © 

A(t) = Aqk sin Qt — — Aqk e cos[(51 — u> prm )t] 

+ ^-AqK€ COs[(il + UJ prm )t] 
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FIG. 2: Visualization of the state space. A no parametric 
excitation. B parametric excitation: Aq = 0.3; e = 0.1; lo = 
0.01755. The lag-time 330, in both cases. 
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FIG. 3: Probability density £ for the distribution of JVo with 
Ao = 0.3; 6 six bound states (levels confined to the well). 
Dashed line - no parametric excitation. Solid line - parametric 
excitation with e = 0.1., parametric frequency u) = 0.01755 
close to the resonant One, UJprm 



The triplet structure could generally have an impor- 
tant bearing on the tunneling in the double-well poten- 
tial. In fact, if a monochromatic pulse at resonant fre- 
quency f2 = Ei — Eq is employed, there is no parametric 
excitation and the dynamics of the occupation probabil- 
ity has the usual form. In contrast, a poor quality non- 
monochromatic pulse may contain the triplet Q ± 0J prm 
so that the deformation of the tunneling dynamics and 
the emergence of the chaotic motion become possible. It 
is worth noting that the amplitude of the driving force, 
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A , can be small so that the frequency shift determined 
by LOprm be tiny. 
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